How to measure the spatial correlations in disordered 
Berezinski-Kosterlitz-Thouless transition ? 
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The Berezinski-Kosterlitz-Thouless transition is a unique two dimensional phase transition, separat- 
ing two phases with exponentially and power-law decaying correlations, respectively. In disordered 
systems, these correlations propagate along lavorable paths, with the transition marking the point 
where global coherence is lost. Here we propose an experimental method to probe locally these 
particular paths in superconducting thin films, which exhibit this transition, and demonstrate theo- 
retically that close to the transition the coherence propagate along a ramified network, reminiscent 
of a percolation transition. We suggest and calculate experimentally accessible quantities that can 
shed light on the spatial correlations in the system as it approaches the critical point. 

PACS numbers: 74.20.-z, 74.78.-w ,74.81.-g 



According to the Mermin- Wagner theorem [l[ there 
can be no long-range order (of continuous symme- 
try) at any finite temperature in two-dimensional sys- 
tems. Nevertheless, one may expect a finite-temperature 
Berezinski-Kosterlitz-Thouless (BKT) transition 
between a phase with power-law decaying correlations 
at T < Tbkt and a phase with exponentially decaying 
correlations at T > Tbkt- This BKT transition is driven 
by the unbinding and proliferation of vortex-anti vortex 
pairs at the critical temperature Tbkt, which destroy 
global phase coherence. It was suggested [f| that such 
a BKT transition can be realized as a superconductor- 
insulator transition (SIT) in thin disordered films, where 
one of the hallmarks of the BKT transition is a universal 
jump in the superfluid stiffness [6| at the critical tem- 
perature. This SIT in thin films has been the subject of 
intense experimental and theoretical research for several 
decades [7[ , yet many puzzles remain. 
In such disordered films, it has been established theo- 
retically [§- 13 1 and observed experimentally 14 - 17 1 that 
the superconducting (SC) order parameter fluctuates 
strongly across the sample, creating "SC islands", where 
the SC order is high, surrounded by areas of weaker 
SC correlations. With increasing temperature the coher- 
ence between neighboring SC islands is suppressed, until 
percolation of coherence from one side of the sample to 
the other is lost, leading to the loss of global SC order 
12L [l3j, [l8[ . This description predicts that local SC order 



may persist even when global SC order is lost, consistent 
with experiments 19l- 24 1. 

It has not been established whether the above process is 
an alternative description of the BKT transition in a dis- 
ordered system or describes a different phase transition. 
The Harris criterion [25| indicates that the BKT tran- 
sition is stable under weak disorder, and recent experi- 
mental |26| and theoretical [13| studies on disordered su- 
perconductors demonstrated the consistency of the tran- 
sition in disordered system with BKT physics (see how- 
ever 12711) . The latter calculation also showed that the 



BKT transition is indeed concomitant with percolation 
of near-neighbor correlations. 

Understanding how coherence propagates in such 
disordered systems is thus crucial to resolving the 
nature of the transition. In this Letter, we propose 
an experimental procedure based on an established 
experimental technique [281430] , whereby the change in a 
global quantity, such as the resistance or the superfluid 
response, is measured as a function of a local perturba- 
tion, such as a charged AFM or STM tip, which can be 
scanned across the system. A significant change in the 
conductance, for example, indicates that the area under 
the tip carries, in the absence of the tip, a significant 
portion of the current through the sample, while small 
changes indicate areas the current tends to avoid. 
Such a procedure produced clear maps of the modes 
of transport in quantum point contacts [2|| |3(| and in 
the quantum Hall regime [29j . Generalizing the same 
procedure to SC thin films will allow spatial resolution 
of the local propagation of coherence in disordered 
superconductors, and, more generally, in disordered 
systems manifesting the BKT transition. (Standard 
scanning STM measurements have already been done 
in thin SC films by eg. [16|, |l7|). Here we produce 



numerical simulations of the experimental procedure on 
a square lattice, and make specific predictions. A main 
result of our work is that indeed certain links carry 
more of the global coherence as the critical temperature 
Tbkt is approached and that these links form a ramified 
network. We then present further evidence supporting 
the percolation hypothesis by showing a connection 
between this network and the Wolff cluster 3l|, which 
is known to percolate at Tbkt 32j. As we demonstrate, 
the measurement will allow direct probe of local and 
long-range correlations in the system, which will provide 
detailed description of the transition. 

Model and Method: A minimal model that cap- 
tures the BKT behavior in two dimensions is the classical 
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XY model. It neglects the fermionic excitations possible 
when cooper pairs are broken into two quasi-particles, yet 
it captures the loss of global phase coherence by a BKT 
transition. This model is relevant to low-dimensional 
superconductors where loss of phase coherence dic- 
tates the transition temperature Tbkt, which is much 
lower than the SC gap, the quasi-particle excitation en- 

ergyiHH. 



The model is described by the Hamiltonian, 



Hxy - - ^2 Ji < 



cos 



(1) 



With Oi G [0, 2ir] an angle in the XY plane, and (ij) 
denotes nearest neighbors. Jij, the random nearest 
neighbor couplings, are taken from a log-box distribu- 
tion, Jij — Joexp(x), with x distributed uniformly be- 
tween —A and +A, and A > thus controls the disorder 
strength. The distribution is chosen such that the me- 
dian value of Jij is Jo, defined as the unit of energy, 
regardless of the disorder A. In what follows we present 
results for A = 5 which provides values of Jij over many 
orders of magnitude while still allowing enough numeri- 
cal accuracy, but we find similar results for a range of A. 
We simulate the model using the classical Monte Carlo 
Wolff algorithm [3l| with periodic boundary conditions 
in both directions. 

The global quantity that we study is the superfluid re- 
sponse, or the helicity modulus. The helicity modulus 
T x in the x direction is the response of the free energy 
F to a phase twist <& x in that direction and is given by 
the expression 



d 2 F 




(2) 



where the sums are over all bonds in x direction and 
N x is the number of such bonds. 0y = 6i — 9j is the 
phase difference between neighboring sites i,j and the 
large brackets denote thermal averages. T y is defined 
similarly for a twist in the y direction, and we define 
T = (T x + T y )/2. At low temperatures, where phase 
fluctuations are suppressed, cos(^) « 1 — \Q~lj-, and the 
XY model is equivalent to a random resistor network. 
Since the resistance of such network with an exponen- 
tial distribution of resistances is given by the percolating 
resistance [IjJ (3(|, which, for the square lattice, is the 
median of the distribution, we expect that T(T) will be 
unaffected by disorder (the value of A) for the log-box 
distributions studied in this work, for low temperatures. 
Indeed we find that T(T) is unmodified by A at very low 
temperatures (T <C min[Jij]) and that Tbkt is reduced 
by less than 10% for large disorder (A = 5). Moreover, 



20 
40 
60 
80 
100 
120 



up 



m 



.■-■*■„'■ ■ '-; 



1 00 1 20 





311118 

20 40 60 80 100 120 



FIG. 1. ST for a typical realization of disorder at A = 5 
for a 128 x 128 periodic lattice, for T = 0.11 -C Tbkt (a) 
and T = 0.89 « Tbkt (b). As temperature increases to- 
wards Tbkt, coherence propagates through specific channels, 
forming a ramified network, (c) The distribution P(ST/T) 
at T — 0.11 — 0.91 showing a constant median (vertical lines) 
and a shift of the weight to the tails with increasing temper- 
atures (change points marked by the arrows). Dashed line: 
the distribution at T — 0.89 divided to the J < 1 (left) and 
J > 1 (right) components (d) P w (probability to belong to a 
Wolff cluster) at T — 0.89 showing the same spatial structure 
of the Wolff clusters as the 8T plot in (b). 



the shape of T(T) is consistent with the BKT behavior, 
i.e. as the system size increases, the (negative) slope in 
T(T) increases for T > Tbkt but not for T < T B kt, 
indicating a finite jump in the large size limit. 

We simulate a local disturbance (such as the STM tip 
in the experiment) by cutting a single bond (i.e. setting a 
single Jij — 0) and then measuring the resulting change 
in T, which we label <5Yy. This procedure is repeated for 
all bonds in the lattice. Removing a bond {i,j} in the 
x— direction results in a change in global helicity modu- 
lus, 



Jij (c 
2pJ K 



<0ij)) + P 4 < sin2 (*«)) 

Jim (sin(%)sin(6>/ m )) 



(3) 



{lm)\\t 



and similarly for a bond in the y— direction. As we are 
interested in linear response, the brackets here denote 
thermal averages with respect to the unperturbed 
Hamiltonian in Eq. [TJ 

Results: Fig. [T] shows a two-dimensional map of <5Ty 
for a typical realization of disorder at A = 5 for a 128 x 128 
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periodic lattice, at two temperatures, T <C Tbkt and 
T « Tbkt- (To display the bond values, we map each 
bond to a site using the established technique of a cover- 
ing lattice [37j ■) Darker bonds mark those bonds that 
cutting them makes a stronger change to the helicity 
modulus. For T <C Tbkt (Fig. [T^) it is evident that 
the bonds that contribute significantly to global coher- 
ence form a dense network. As T approaches Tbkt (Fig. 
[TJb) , these bonds form a very ramified network, consisting 
of distinct channels, even though the disorder does not 
have any spatial correlations. In Fig. [1}; we plot the spa- 
tial histogram of the values of ST(T)/T(T) for a range of 
temperatures. Three features in the histogram are note- 
worthy: (i) the median of <5T/Y does not change with 
temperature; (ii) As temperature is increased towards 
Tbkt, the peak of the distribution (the range of ST /T 
between the arrows) gets suppressed and the weight is 
shifted to the tails of very strong and very weak ST val- 
ues, (iii) The distribution is a sum of two peaked dis- 
tributions, one centered around high ST (near the right 
arrow) and one around small values (left arrow). The two 
broken lines denotes the distributions of ST/T generated 
by the bonds with a microscopic Jij < 1 (left) and those 
with > 1 (right), revealing a correlation between ST 
and (see below). As disorder (A) is increased, the 
peaks grow farther apart. 

Consequently, the following physical picture emerges: at 
low temperature the sample is mostly coherent and there 
are many routes for the global phase coherence to propa- 
gate through the sample, thus many bonds give a signif- 
icant contribution when cut. As temperature increases 
towards Tbkt, larger regions lose coherence, so the co- 
herent routes gradually thin out, limiting the propagation 
of global coherence to a smaller set of bonds, forming a 
ramified network of channels. 

Naively, one might suspect that there exists a trivial re- 
lation between the bare microscopic couplings Jij and 
their global effect STij, as might be inferred from Eq. 
[3l Yet Fig. [2^i reveals a more interesting relation: for 
sub- median bonds Jy < 1, there is a linear correlation 
ST^ Jij as suggested by the unit slope of the log 
plot. However, for half the bonds of the lattice, for which 
Jij > 1, STij is only weakly dependent on J,j (the exact 
relation depends on A). The reason is that these bonds 
belong to the percolating coherence cluster, and the ef- 
fect of cutting such a bond on ST^ depends no longer 
only on its bare coupling Jij but rather its location in 
the network. We consider this behavior and its onset 
at the median bond (which is the percolation threshold 
for the square lattice) as further evidence supporting the 
percolation scenario. 

In order to shed more light upon the spatial structure of 
the correlation, we turn our attention to the Wolff clus- 
ters 31(. A priori, these are mere numerical constructs 
that were introduced to overcome the numerical problem 
within the Metropolis Monte Carlo algorithm, which is 



based on updates of local variables (such as spins), de- 
pending on Boltzman weights. Due to the divergence of 
length scales near the transition, spins are highly cor- 
related on large scales, and single spin fli ps a re highly 
improbable. To address this problem Wolff [3l| proposed 
a very efficient algorithm which eliminated this critical 
slowing down by allowing clusters of spins, called the 
Wolff clusters, to flip together. These clusters are con- 
structed stochastically according to near-neighbor corre- 
lations. It was later suggested [38| and then proven rig- 
orously [32| that the efficiency of the algorithm lies in its 
correct capture of the physics of the BKT transition. The 
size of the Wolff clusters, that reflect the correlated spins, 
grows towards the BKT transition, the latter happening 
exactly where the Wolff cluster percolates through the 
whole system. Note that this proof applies for a clean 
system - the random structure of the Wolff cluster is a 
result of the stochastic nature of its construction, so one 
can imagine that they reflect the random temporal corre- 
lations of the spins. Nevertheless, due to the uniformity 
of the system, the probability of every bond < ij > to 
belong to a Wolff cluster, Pjj is the same. This situation 
is very different in disordered systems. The Wolff clus- 
ters, generated during the Monte Carlo calculation, are 
highly anisotropic and reflect the spatial disorder in the 
system. In Fig. QJd) we display a two dimensional plot 
of the probability P^j at T w Tbkt- Indeed, P# and ST 
(Fig. [Tfb)) exhibit very similar spatial structure, so the 
bonds with high ST, darkest in Fig. [TJb), are also the 
ones with the highest probability to be part of a Wolff 
cluster, brightest in Fig. [Tf d) . Hence the bonds that 
support the propagation of coherence (i.e have a large 
ST) are also, as expected, the ones that have long range 
correlations, as reflected by the Wolff clusters. This ob- 
servation is further supported by FigJSJb), which displays 
a normalized scatter plot of ST^ vs P^j, demonstrating 
the clear correlation between over ten orders of magni- 
tude. If the proof that the Wolff clusters percolate at 
the BKT transition extends to the disordered case, as is 
indicated numerically, it is yet another support for per- 
colating nature of coherence at the BKT transition. 
A further check on the changing spatial structure with 
temperature can be revealed by studying the correlation 
function 



Cr(R) 



E 



<ij>-<kl> 



=R (ST^ - ST) (6T kl - ST) 



E 



\<ij>— <kl>\=R 



(4) 

where the sum is over all pairs of bonds, distance R apart, 
and ST is the asymptotic average of STij, so that CV(i?) 
decays to zero at large R. A similar correlation func- 
tion Cw(R) can be defined in terms of Py. Due to 
the high correlations between P-j and STij (FigEJb)) 
we find a similar behavior for these two correlation func- 
tions, though Cw(R) displays smaller statistical fluctua- 
tions. FigHIc) depicts Cw{R) at a range of temperatures 
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FIG. 2. (a) Normalized scatter plot of In ST vs. In J show- 
ing correlation only for the sub-median bonds (In J < 0). (b) 
Normalized scatter plot of ln<5Y vs. lnP^ (probability to be- 
long to a Wolff cluster) showing a clear linear correlation, (c) 
Cw(R), the correlation function for P w for a range of temper- 
atures on a linear and (d) log-log scale. The correlation decays 
more slowly at higher temperatures, indicating the sparsity of 
the coherence propagation. 



and Figl^d) depicts it on a log-log scale. At low tem- 
peratures, as the network supporting coherence is dense, 
Cw{R) decays quickly to zero. As temperature increases 
there is a marked slower decay of Cw (R) , demonstrating 
the ramified nature of that network. For the system sizes 
we simulated, the correlation function doesn't follow any 
simple decay law, possibly due to finite size effects. 

Summary and discussion: The effect of disorder on 
the BKT transition is a hotly debated problem in con- 
densed matter physics, and its relation to percolation, 
or the emergence of a new critical point are open is- 
sues. In this work we propose an experiment to probe 
the spatial structure of global phase coherence in two di- 
mensional superconductors, that are expected to display 
such a transition. Our proposition is based on similar 
experiments that probe conductance channels through 
quantum point contacts and in quantum hall systems 
(while, for theoretical convenience, we display results for 
the phase stiffness, we expect similar results for the cur- 
rent). We present numerical results for a minimal model 
that (in its clean version) displays a BKT transition, in 
the strong disorder regime. Our numerics point at an 
intriguing mechanism where the global coherence passes 
predominantly through specific channels that can be ex- 
perimentally probed. The correlation functions we have 
calculated, which can also be measured experimentally, 
also support the idea that as the critical temperature is 
approached, the coherence network becomes more rami- 
fied, with percolation lost at Tbkt- This picture is fur- 



ther supported theoretically by the robust correlation be- 
tween these channels and the structure of the Wolff clus- 
ters, objects known to percolate at the phase coherence 
transition of clean systems. Moreover, our tests at A > 5 
show that the two peaks (dashed lines in Figfjjc) move 
apart as disorder increases. This suggests the tempting- 
possibility that at large enough A the two peaks separate 
out and can thus give the experimentalist a clear map 
of where coherence propagates, in addition to a direct 
determination of which of the local Josephson coupling 
are below the median as opposed to above the median. 
It is interesting to note that as temperature increases to- 
wards Tbkt, a significant fraction (up to approx. 40%) 
of the ST values become negative, indicating an increase 
in coherence when these bonds are cut. Such behavior is 
certainly possible in coherent systems, but is a departure 
from known results on random resistor network 3^, 40j ]. 
which is relevant to the XY-mode\ at low temperatures. 
Again, this observation is consistent with the percolation 
pictures: cutting links that do not participate in the net- 
work that support the propagation of coherence, not only 
does not impede the global coherence, but may, in fact, 
increase it by giving more weight to the percolating net- 
work. The experiment suggested here and the numerical 
method can, in fact, be applied to any physical system 
involving a phase transition driven by phase fluctuations, 
and may also be relevant to high-Tp superconductors in 
part of their phase diagram. [331 ] . 
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